Competition between the BCS superconductivity and ferromagnetic spin fluctuations 

in MgCNi 3 
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The low temperature specific heat of the superconductor MgCNi3 and a non-superconductor 
MgCo.ssNis is investigated in detail. An additional contribution is observed from the data of MgCNi3 
but absent in MgCo.ssNis, which is demonstrated to be insensitive to the applied magnetic field even 
up to 12 Tesla. A detailed discussion on its origin is then presented. By subtracting this additional 
contribution, the zero field specific heat of MgCNi3 can be well described by the BCS theory with 
the gap ratio (A/fcsT c ) determined by the previous tunneling measurements. The conventional 
s-wave pairing state is further proved by the magnetic field dependence of the specific heat at low 
temperatures and the behavior of the upper critical field. 
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Since the discovery of the new intermetalic perovskite 
superconductor MgCNi 3 pj, plenty of efforts have been 
focused on the superconducting pairing symmetry in 
this material because its conduction electrons are de- 
rivedpredominantly from Ni which is itself a ferromag- 
net 0, 0, S @- However, up to now, there is still 
not a consensus on this issue. The measured penetra- 
tion depth 0, critical current behavior Q and earlier 
tunneling spectra Q suggested an unconventional super- 
conductivity, the later tunneling data supported the 
s-wave pairing symmetry and gave a reasonable inter- 
pretation on the contradiction to the result in Ref 8]. 
The s-wave pairing has also been demonstrated by the 
13 C NMR experiments [l0| and the specific heat mea- 
surements HHEl El- To our knowledge, 
all the previous reports on the specific heat of MgCNi 3 
ED- lJ HI LUi IR* UjJ were characterized in the framework 
of a conventional phonon-mediated pairing. However, 
there is an obvious deviation of the experimental data 
from the prediction of BCS theory in the low tempera- 
ture i- e -> the entropy conservation rule is not satis- 
fied. Such deviation has been interpreted by the presence 
of unreacted Ni impurities in Refsj^, H^ - whereas it is 
still prominent in the samples without Ni impurities |14| . 
On the other hand, strong spin fluctuations have been 
observed in MgCNi3 by NMR experiment [HJ, which is 
suggested to be able to severely affect the superconduc- 
tivity in MgCNi 3 @, H liT lfTll [Icj or even induce some 
exotic paring mechanism [2J. Consequently, the behav- 
ior of the specific heat will inevitably be changed by the 
spin fluctuations. Therefore, before a real pairing mech- 
anism being concluded from the specific heat data, we 
have to carefully investigate how the ferromagnetic spin 
fluctuations contribute to the specific heat of MgCNi3 . 

In this work, we elaborate on the specific heat (C) 
of MgC^Nia system both in normal state and supercon- 
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ducting state. A low temperature upturn is clearly dis- 
tinguished in the C/T vs T 2 curves and found to be in- 
sensitive to the applied magnetic field. By doing some 
quantitative analysis, we present the evidence of most 
possible mechanisms responsible for this upturn. Af- 
ter subtracting this additional contribution, a well de- 
fined BCS-type electronic specific heat is extracted. The 
temperature dependence of the upper critical field and 
the field dependence of the low temperature specific heat 
also supports such conventional BCS superconductivity 
in MgCNi3. These analyses indicate that although the 
spin fluctuations may suppress the pairing strength in 
MgCNi 3 , the superconductivity is certainly not induced 
by any exotic mechanism. 

Poly-crystalline samples of MgC^Nia were prepared 
by powder metallurgy method. Details of the prepara- 
tion were published previously |17| . The superconductor 
MgCNi3 has a T c of Q.IK and the non-superconductor 
MgCo.ssNis was synthesized by continually reducing the 
carbon component until the diamagnetism was com- 
pletely suppressed. The heat capacity data presented 
here were taken with the relaxation method 01 based on 
an Oxford cryogenic system Maglab in which the mag- 
netic field can be achieved up to 12 Tesla. Details of the 
sample information and the measurements can be found 
in recent report It should be emphasized here that 
the Cernox thermometer used for calorimetry has been 
calibrated at 0, 1, 2, 4, 8 and 12 Tesla, and the calibration 
for the intermediate fields is performed by an interpola- 
tion using the result of the adjacent fields. Therefore, any 
prominent field dependence of the specific data should re- 
flect the intrinsic properties of the measured sample. 

In general, the low temperature specific heat C(T, H) 
of a superconductor consists four main contribution s by 
neg lecting the component of the nuclear moments [2(1 
|29| . each has a different dependence on T and two of 
which depend on H , also in different ways, 

C(H,T) = C mag (H,T) + C DOS {H, T)+ 1o T + C ph (T) 

(1) 
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FIG. 1: The low-temperature specific heat of MgCNi3 at 
various magnetic fields from OT to 1 2T. The thick solid line 
denotes the universal background including all the normal 
state data for different fields. 



where 70T represents a spare zero-field linear term and 
C p h is due to the lattice or phonon contribution. The 
Debye phonon specific heat C p h — (3T 3 can usually de- 
scribe the lattice contribution at low temperatures. How- 
ever, the departure from T 3 -law has often been observed, 
which is due to the fact that the density of modes of 
the phonon in real solid does not follow the assumed 
w 2 -law, here u> is the angular frequency of a harmonic 
wave associated with the lattice vibration. In such case, 
the deviations may be expanded in higher order terms 
such as T 5 , T 7 , etc.. The ii-dependent terms in Eq. Q), 
i.e., C mag (H) and Cdos(H), are the contributions as- 
sociated with magnetism and the electronic density of 
states ( DOS ), respectively. If there is no magnetism 
associated contribution, the normal state specific heat 
at low temperature can be approximatively described as 
C n (T) = 7„T -(- /3T 3 in the framework of metal theory. 
Therefore, a linear relation can be obtained by plotting 
the normal state data as C n (T)/T vs T 2 , and its inter- 
cept and slope correspond to j n and [3, respectively. 

The low temperature specific heat at various magnetic 
fields up to 12 Tesla is plotted as C{T)/T vs T 2 in 
Fig. ^ Two important features should be emphasized 
here. First, all the normal state data at various mag- 
netic fields merge into one which is consistent with 
the results reported by other groups 0, 0. Second, 
this common normal state background remarkably devi- 
ates from the linear relation as discussed above. In Fig. [21 
only the normal state data are re-plotted in an magnified 
scale. In order to survey the normal behavior at very 
low temperature, the magnetic field up to 12 Tesla was 
applied in measurements, which exceeds the highest up- 
per critical field of our sam ple and is 4 Tesla higher than 
that used by other groups [lj, [l5| . The specific heat of 
the non-superconductor MgCo.ssNis is also presented in 
Fig. |21 as a comparison. It is obvious that the perfect lin- 
ear relation of C(T)/T vs T 2 is satisfied for MgC .85Ni 3 , 
which is a striking contrast with the case of MgCNi3. 
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FIG. 2: A plot of C/T vs T 2 for MgCNi 3 and MgCo.ssNis in 
different magnetic fields. The data of MgCo.8sNi3 can be well 
defined by a straight line while that of MgCNi3 remarkably 
deviate the linearity. The dashed line is a linear extrapola- 
tion of the high temperature data of MgCNi3 and the upper 
solid line is the theoretical fit considering higher order phonon 
contributions. 



The obvious upturn in the low temperature C(T)/T 
vs T 2 curves of MgCNi3 can not be associated with Ni 
impurites since the X-ray diffraction pattern shows no 
indication for Ni impurites [ll|. To say the least, if there 
is still extreme small content of Ni impurites leading to 
the prominent low temperature upturn of C/T, the field 
dependence of its specific heat should also be obvious, 
which is clearly inconsistent with our experimental re- 
sults. Moreover, if this upturn is due to the excess free 
Ni in MgCNi3, it should also be observed in MgCo.ssNis 
because of the similar process of synthesizing these two 
samples. Quantitatively, taking the data from references 
|2lLl22| yields for 10% of superfluous Ni an upturn which 
is at least two orders of magnitude smaller than the ob- 
served one. Therefore, the contribution of the excessive 
Ni can be neglected comparing with the whole specific 
heat. Furthermore, the possible Schottky anomaly is 
presented in Fig|31 its field dependence is obviously too 
strong to compare with our experimental result ( nearly 
field independent ). 

It is found that this upturn can be well fitted if the 
above mentioned T 5 term is considered ( see the up- 
per solid line in Fig. El )• In other words, the depar- 
ture from the T 3 behavior may be due to the non-Debye 
phonon DOS, which is consistent with the notable dif- 
ference of the Debye temperature between MgCNi3 and 
MgCo.ssNis 01 ■ If t ne electron-phonon coupling is in- 
deed the origin of superconductivity in MgCNi3, it is 
reasonable to associate the disappearance of supercon- 
ductivity in MgCo.ssNis with the remarkable difference 
of its phonon DOS from that of MgCNi3. However, some 
careful work is needed to understand such obvious dif- 
ference of phonon structure between these two samples, 
since they have similar crystal lattices and chemical com- 
ponents. 
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FIG. 3: Comparison between the field-insensitive low- 
temperature upturn in the specific heat of MgCNi3 ( AC = 
C — 7T — (3T 3 ) and the calculated field dependent Schottky 
anomaly. 



FIG. 4: Electronic specific heat 7„(T) versus T of MgCNi 3 in 
the normal state ( three different f3- values are selected in order 
to avoid artificial errors. ). The solid lines are theoretical fits 
to spin-fluctuation model. All the curves except the top two 
are shifted downwards for clarity. 



Another possible explanation of the above mentioned 
low temperature upturn is the existence of strong spin 
fluctuations due to the higher DOS at fcrmi energy ( 
N(E F ) ) of MgCNi 3 than that of MgC . 85 Ni 3 [lj, con- 
sequently, the coupling between the electrons and spin 
fluctuations in MgCNi3 should also be stronger. The fer- 
romagnetic spin fluctuations have been demonstrate d by 
NMR experiments [13. Doniach and Engelshberg [23] 
and Berk and Schrieffer [24| showed that the absorp- 
tion and re-emission of spin fluctruations renormalizcs 
the electronic self-energy, leading to an enhanced effec- 
tive mass at low temperatures. This effect manifests it- 
self as a low-temperature enhancement of the electronic 
specific-heat coefficient, X s f, which depends on tempera- 
ture as T 2 ln(T/T s f) (here T s f is the characteristic spin- 
fluctuation temperature) at low temperature. Consider- 
ing the presence of ferromagnetic spin fluctuations, the 
normal state specific heat of MgC^Nia can be expressed 
as follows, 

C n (H = 0,T) = A[l + X ph + X sf (T)}T + l0 T + l3T 3 (2) 

where /3T 3 are the contributions of phonon excitations, 
X s fT and X p hT represent the contributions of effective 
mass renormalization due to the electron-spin fluctua- 
tion coupling and the electron-phonon coupling, respec- 
tively, and A is a constant correlated with N(Ep). It 
can be seen from Eq. that the deviation from the 
linear dependence of C{T)/T on T 2 is due to the tem- 
perature dependence of X s f. Moreover, Beal-Monod, Ma, 
and Fredkin [2f| have estimated the shift 8C/T caused 
by an applied field H to be 



0.1 ( 



Mg )2 S 
ksTsf InS 



(3) 



where S is Stoner factor. Eqs. J5J and © indicate 
that the possible magnetic field dependence of the nor- 
mal state specific heat is completely determined by the 



spin fluctuations. For simplicity, Eq. (J2J can be rewrit- 
ten as C n (H = 0, T) = 7„(T)T + /3T 3 , in which j n (T) = 
A[l + X ph + X sf (T) + 70 /A}. Therefore, the 7„ ~ T rela- 
tion directly reflects the temperature dependence of A s / . 
In Fig. 0] we present the determined 7 n (T) by selecting 
various (3- values. Fitting the 7n(T) relations to the for- 
mula of A(l + BT 2 ln(T /T sf )) yields T sf varying from 
13 to 16K. By inserting the determined T s f, calculated 
Stoner factor S 2] 

and the highest field value in our mea- 
surements into Eq. J3J), we can estimate the shift SC/T 
caused by the applied field to be less than 2%, which is 
in agreement with our experimental results. However, 
if this explanation is correct, we must understand the 
collapse of the entropy conservation around T c caused 
by considering such additional electronic specific heat, as 
discussed below. Therefore, the specific-heat contribu- 
tion of the spin fluctuations themselves may be another 
candidate responsible for the low temperature upturn in 
specific heat of MgCNi3 . 

Despite the true mechanism of the low-temperature 
upturn of C /T, this additional specific heat contribution 
should be regarded as a part of the normal-state back- 
ground of the superconducting specific heat below the up- 
per critical field H r o (T). In earlier analysis to the specific 
heat data 0, ITTL Il2| , this additional part of background 
has been neglected more or less below H C 2(T). We point 
out here that neglecting this additional contribution will 
lead to the collapse of the entropy conservation as re- 
ported in references 0, Il5j| . This opinion is motivated 
by the subsequent analysis. As shown in Fig. 0a), the 
normal state background ( as shown in Fig. [21 ) has been 
subtracted from the zero-field specific heat data, the en- 
tropy difference AS(T) = J Q T dT(AC/T) is presented in 
the inset of Fig. 0Ja), here AC = Ch=o — C n . It is 
found that the entropy conservation is then well satis- 
fied, indicating that the remainder is the contribution of 
superconducting state. Such analysis has also been ap- 
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TABLE I: Fits to BCS model for zero-field specific heat. 
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FIG. 5: Fitting the specific heat data (Ch=o — C„) mea- 
sured by different groups [U Q to BCS model. The 
experimental data are denoted by open circles and the fits to 
a-model are denoted by solid lines. The dotted lines are fits 
to the original BCS model. The inset in (a) is the entropy 
difference by integration of AC/T as presented in (a). 



plied on the data measured by Walte et al. 0] and Lin 
et al. 0| , respectively, the entropy is also conserved and 
the low temperature anomaly as mentioned in Ref.|l5| 
completely vanishes. 

When the superconductivity in MgCNi3 is investi- 
gated, the spin fluctuations can not be neglected because 
it may compete with superconductivity j^, 1 1 1| or even 
lead to an exotic pairing mechanism other than the con- 
ventional s-wave |2j ■ If the phonon-mediated pairing wins 
in competing with the spin fluctuations and hence the 
effect of the spin fluctuations only suppress the electron- 
phonon coupling; ( or pairing-strength ) [ijj , the so-called 
a-model |2g based on BCS theory should be a good 
choice to describe the measured thermodynamic param- 
eters. Comparing with the original BCS-model, the only 
adjustable parameter in this a-model is the gap ratio 
A(0)/fcsT c . This model has been successfully applied 
to strong coupling systems such as Pb and Hg. Then 
we try to fit the superconducting part of the zero-field 
specific heat to the a-model, the results are presented 
in Fig. 03 All the data can be well described by this 
revised BCS model with the best fitting parameters ( 
i.e., A and T c ) listed in Table. [I] These fits yield a 
gap ratio A(0)//cbT c w 2.06, corresponding to the max- 
imum gap A(0) « 1.2meV which is in good agreement 
with our previous tuneling measurements 0. From the 
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above discussions, we can conclude that the coexistence 
and competition of spin fluctuations and phonons does 
not change the phonon-mediated pairing mechanism of 
MgCNi 3 . 

In order to further verify this picture, we investigate 
the field dependence of the low temperature specific heat 
of MgCNi3. It is known that the electronic specific 
heat in magnetic fields can be expressed by C e ; (T, H) — 
C ei {T, H = 0) + j(H)T. The magnetic field dependence 
of j(H) is associated with the form of the gap function 
of the superconductor. For example, in a superconduc- 
tor with line nodes in the gap function, the quasiparticle 
DOS ( N(E) ) rises linearly with energy at the Fermi 
level in zero field, N(E) cx \E — Ep\, which results in 
a contribution to the specific heat Cdos — olT 2 |27j . 
In the mixed state with the field higher than a certain 
value, the DOS near the Fermi surface becomes finite, 
therefore the quadratic term Cdos = a T 2 will disap- 
pear and be substituted by the excitations from both 
inside the vortex core and the de-localized excitations 
outside the core. For d-wave superconductors with line 
nodes in the gap function, Volovik et al. pointed out 
that in the mixed state, supercurrents around a vortex 
core cause a Doppler shift of the quasi-particle excita- 
tion spectrum. This shift has important effects upon 
the low energy excitation around the nodes, where its 
value is comparable to the width of the superconducting 
gap. For H » H cl , it is predicted that N(E F ) oc H 1/2 
and C D os = &y(H)T = ATH 1 / 2 at low tempera- 
tures |28|]. This prediction has been well proved for 
hole doped cuprates [2jJ. Whereas in a conventional 
s-wave superconductor, the specific heat in the vortex 
state is dominated by the contribution from the local- 
ized quasi-particles in the vortex cores. From the Bo- 
goliubov equations assuming noninteracting vortices, the 
DOS associated with the bound excitations is derived as 
N(E) oc B(H) |3(j hence the contribution of the vortex 
cores to the specific heat is Cdos °c B{H)T 0,H2- It 
is also theoretically derived that the experimentally ob- 
served downward curving C(H) is caused by the flux line 
interactions near H c \ and the possible expansion of the 
vortex cores 0, . According to the above discussions, 
the specific heat coefficient j(H) of conventional s-wave 
superconductor should linearly depend on the magnetic 
field well above H c \. 

Fig. shows the field dependence of j(H) — 7(0) of 
MgCNiq below 3K. The data reported by different groups 
[ill [L4 fl5| merge into each other by timing a prefactor 
A close to unity. It is found that 7 linearly depends on 
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FIG. 6: The magnetic field dependence of the specific heat 
coefficient ^f(H) — 7(0) at low temperatures. 
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FIG. 7: The comparison of the temperature dependence of 
the upper critical field H C 2(T) with the BCS-like descriptions 
in lower and higher temperature limits. 



H above 0.5T and persists up to 8T which is close to 
the upper critical field of MgCNi3. The legible linear- 
ity of A7 ~ H relation at higher field and its negative 
curvature below 0.5T are in good agreement with the 
above mentioned behaviors of conventional s-wave super- 
conductors. It may be argued that the low temperature 
limit of about 2K in our measurements is not low enough 
to distinguish the ci-wave's A7 ~ iJ 1 / 2 -law. However, it 
should be emphasized that the observed A7 ~ H rela- 
tion is nearly universal at low temperatures below upper 
critical field, which is very similar to the behavior of V^Si 



[33, a typical conventional s-wave superconductor. 

Finally, we compare the temperature dependence of 
the upper critical field H C 2(T) with the prediction of BCS 
theory in which the H c2 (T) can be expressed as follows, 



H C 2{T) w 1.74iJ c (0)(l - T/T c ) (T c - T < T c ) (4a) 



H c2 {T) w ff c2 (0)[l - 1.06(T/T C ) 2 ] {At low T) (4b) 
As shown in Fig. [7| the best fitting to BCS model is 
denoted by solid lines. At lower temperature, the ex- 
perimental data can be well described by Eq.(0J)). For 
the higher temperature near T c , a prefactor of 1.65 is 
obtained instead of the theoretical prediction of 1.74 as 
expressed in Eq.l^). Nonetheless, the BCS model is still 
a preferred description for H c2 (T) of MgCNi 3 considering 
the stronger electron-phonon coupling and the presence 
of ferromagnetic spin fluctuations. 

In summary, we have investigated the specific heat 
data of MgC^Nia system. A remarkable field independent 
contribution is found in MgCNi3, reflecting the departure 
of normal-state specific heat from T 3 -law. By removing 
this contribution, the zero-field data is well described by 
the a-model ( a slightly revised BCS model ). The con- 
ventional s-wave superconductivity is further supported 
by the linear field dependence of specific heat coefficient 
7(H) and the BCS-like temperature dependence of up- 
per critical field H C 2(T). It is then concluded that, al- 
though electron-magnon ( spin fluctuations ) coupling 
coexists and competes with electron-phonon coupling ef- 
fect in MgCNi3, it only acts as pair breakers while does 
not induce a new exotic superconductivity. 

Note added: Most recently, the carbon isotope effect in 
superconducting MgCNi3 observed by T. Klimczuk and 
R.J. Cava indicates that carbon-based phonons play a 
critical role in the presence of superconductivity in this 
compound 33]. 
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